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Abstract 

It is shown that thermal fluctuations present in a simple non-degenerate relativistic fluid 
satisfy a wave equation in the Euler regime. The characteristic propagation speeds are cal- 
culated and the classical expression for the speed of sound is recovered at the non-relativistic 
limit. Implications and generalizations of this work are analyzed. 

1 Introduction 

Relativistic transport theory has dramatically increased its interest due to the detection of high 
temperature plasmas produced in the Relativistic Heavy Ion Collider (RHIC). In this context, 
an Euler fluid description provides a good approximation for events involving Au-Au collisions 
[T]. Generalizations involving dissipative effects have been developed taking into account recent 
experimental data The first works regarding relativistic hydrodynamics can be tracked down 
to the pioneering 1940 Eckart's monographs [3], and to the relativistic fiuids section included in 
the Landau-Lifshitz fiuid mechanics textbook [4]. The explicit form of the linearized transport 
equations obtained within Eckart's framework raised serious doubts concerning the stability and 
causality properties of the system [5] [6]. Indeed, it was only recently observed that the so-called 
stability problem in relativistic hydrodynamics is due to the heat-acceleration coupHng introduced 
in Eckart's work [7J. Following these ideas, it became pertinent to revise the causality properties 
in this kind of systems, focusing in the possibility of generating a hyperbolic partial differential 
equation describing temperature fluctuations. This work tackles the problem for an Euler fluid and 
suggests some new insights for this issue while examining the linearized equations in the Navier- 
Stokes regime without resorting to extended formalisms. 

In section 2 the basic formaHsm is presented on the basis of relativistic kinetic theory for an 
inert dilute fluid, emphasizing the role of the Enskog transport equation. Section 3 is devoted to 
the analysis of the linearized transport equations in the Euler regime, by means of a derivation 
of a wave equation describing thermal fluctuations in the relativistic case, and its corresponding 
non-relativistic limit. Some flnal thoughts regarding the causal properties of relativistic fluids in 
the dissipative case are included in the final section of this work. 
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2 Kinetic foundations and transport equations 



For decades, relativistic kinetic theory has been successfully applied in the study of high temperature 
fluids jl]. The starting point here is the relativistic Boltzmann equation for a simple fluid in the 
absence of external forces: 



VV,C. = Jiff) (1) 

In Eq. idl), / is the distribution function in the phase space, J{f f) is the collisional kernel, and 
the molecular four velocity, is given by 

v"' = {rfw\"fc) (2) 
where w is the molecular velocity (three spatial components). As usual, 7=11 — ^^^-^ 1 . All 



latin indices run from 1 to 3 and the greek ones run up to 4. A signature (1, 1, 1, —1) is taken, so 
that u°'Ua = — c^. The relativistic generalization of Enskog's transport equation can be casted in 
the form |9]: 

^(n(vI/)) + (n(u/*)),-0 (3) 
where n is the particle number density and the average of the collisional invariant (vf) is defined as 

{^) = ^Jl^fdv* (4) 

with dv* = 7^^^ [ig. 

In the Euler regime, all averages are calculated using the equilibrium (Juttner) distribution 
function, valid for a non-degenerate gas [8]: 

where = {v^) is the hydrodynamic velocity, z = is the relativistic parameter and K2 is 
the modified Bessel function of the second kind. Derivatives with respect to can be explicitly 
evaluated in Eq.|[5]). After this operation, for the sake of simpHcity, all calculations will be performed 
in the comoving frame of the fiuid. 

Now, the collisional invariants are 4" = 1 (a constant), mw^j (the three-momentum) and mc^7 
(the mechanical energy). For "ii — 1 the continuity equation follows immediately 

^ (n) + {nu% = (6) 

Substituting ^ = mw'''y, the momentum balance equation is obtained: 

d 

— (n (mw^-f)) + {nm {vJ-w^-f)),^ = (7) 
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The use of Eqs. ^ and ([5]) allows to rewrite Eq. ^ in terms of the local thermodynamic variables: 



-(.s+,)-(u')+fcn^ + ^T^ = (8) 



c 

where the internal energy per particle e reads: 

e ^ 2,nkT + nmc^ (9) 
and the pressure satisfies the sate equation 

p = nkT (10) 
Finally, for ^ = mc^j, the resulting balance equation reads: 

^ (n(7mc2)) + {n{w\jmc'))),^^0 (11) 



or, in terms of the thermodynamic variables: 

d{ne) 
dt 



pe = (12) 



in Eq. (fT2l) we have defined 6 — u"^ . The set of equations i|6l8ll2p is highly nonlinear and its full 
treatment is rather complex. For a system close to equilibrium we shall Hnearize this set in order 
to perform a fiuctuation analysis for the thermodynamical variables. 

3 Linearized equations and causality analysis 

In order to proceed with the analysis of the Euler system I|6l8ll2p close to equihbrium, we decompose 
any thermodynamical variable X into a constant average value Xq and a space and time dependent 
fiuctuation SX, so that 

X = Xo + SX (13) 

According to this definition, neglecting second order terms, the Hnearized continuity equation, 
obtained from ^ reads: 

^^{6n)+noSe^0 (14) 
Analogously, the Hnearized momentum balance for the longitudinal mode 59 becomes: 

Po^^iSe) + UokV'^iST) + kToV'^iSn) = (15) 

where we have defined po = "°^°+p° . For the linearized energy balance equation we get: 
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Here, the heat capacity (per particle) is given by 

—) (17) 

In order to decouple the system and establish a partial differential equation for ST, we first solve 
for SO in both sides of Eqs. lfT4|) and lfT6|) . Equating the results we obtain the useful relation: 

We derive with respect to time in both sides of Eq. lfT5|) : 

Pog^^^ + nokV^g-^iST) + kT,V^-{5n) = (19) 
so that, inserting the expression for ^Sn from equation (fT8|) we obtain: 

P.^SO + n.kV^-iST) + knV^ -iST) ^ (20) 

One first time derivative is immaterial in each term of Eq. ((20l) . so that after an elementary 
arrangement of terms we can write the wave equation for thermal fiuctuations: 

Thus, the propagation speed (C^) of a thermal wave in a relativistic Euler fiuid is: 

q2 _ npToicv + k)k ^22^ 

Po^v 

In the non-relativistic limit, z ^ 0, Cy |A: and po Uom, so that (|22|) yields 

- l^c' (23) 

which is the non-relativistic speed of sound. Also, the relativistic propagation speed l(22|) can be 
rewritten in terms of z as: 

^2 _ + k)z n , . 

Cy y^Z -h y 

It is interesting to notice that some authors perform a similar analysis for Sn neglecting tem- 
perature fiuctuations, and only taking into account Eqs. ([6]) and ([H, in order to establish a wave 
equation for density fiuctuations [H]. In that case it is immediate to find out that the corresponding 
propagation speed is: 
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Figure 1: Comparison of fluctuation propagation speeds for the full relativistic case(solid), 

non-relativistic case (long dashed), Sn fluctuations neglecting thermal fluctuations 

(short dashed) and ST fluctuations neglecting number density fluctuations (dotted). 
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In the same order of ideas, one can make a simple analysis neglecting the number density 
fluctuations and taking into account only Eqs. ([8]) and lfT6|) . In this case, the expression for a wave 
equation for thermal fluctuations reads: 

6-2 = ^ . , .X (26) 

k K2(l/z)) 

Figure 1 shows a comparison of the characteristic speeds for increasing z. 



4 Final remarks 

It has recently been proved the nonexistence of generic instabilities in the linearized transport 
equations at the Navier-Stokes regime [9]. In this paper it is shown that, in the Euler regime, there 
is no causality problem. The linearized transport equations become a hyperbolic system and, for 
further research, it can be taken as a starting point for a simplified calculation and for vaHdation 
of numerical work in the non-linear case. 

The non-relativistic Hmit has been recovered, as expected, and thermal fiuctuations also satisfy 
a hyperbolic partial differential equation. In most textbooks, the establishment of the (parabolic) 
heat equation is based on an extension of Eq. (fT2l) including heat conduction, neglecting velocity 
fluctuations. On the other hand, if the linearized equation of motion (fT5|) is taken as the basis of the 
description of thermal fluctuations, then a causal equation is obtained for the non-dissipative fluid. 
Thus, for the dissipative case it is suggested that the suitable generaHzation of the whole Hnearized 
system l|6|8|12p should be taken into account, emphasizing the role of Eq.|[8l) when analyzing causal 
properties of the system. Neglecting velocity fluctuations clearly leads to non- causality. It can 
be noted, also, that density fluctuations (neglecting the thermal ones) and thermal fluctuations 
(neglecting the density ones) present different propagation speeds, satisfying the relation C^ + C^ = 

. Moreover, taking S9 = is unrealistic, since when the fluid is at rest the mean velocity is zero. 
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but the fluctuations do not vanish. The approximate expressions (|25|) and ((26l) may be useful in 
particular situations involving dissipative effects. This opens a line for future research. 

The authors wish to thank A.L. Garcia-Perciante for her valuable comments for this work. 



References 

[I] See for example P.P. Kolb and U. Heinz Hydrodynamic description of ultrarelativistic heavy-ion 
collisions in Quark Gluon Plasma 3, R.C. Hwa and X. and N. Wang Eds., World Scientific, 
Singapore pp 634-714 (2004) |arXiv:nucr-th76305684] . 

[2] A. K. Chaudhuri, Phys Rev C 74, 044904 (2006)[ ajXiv:nucl-th/06040T^ . 

[3] C. Eckart, Phys. Rev. 58, 267 (1940); ibid 58, 919 (1940) 

[4] L. Landau and E.M. Lifshitz, Fluid Mechanics Addison Wesley, Reading Mass. (1958). 
[5] W. A. Hiscock and L. Lindblom; Phys. Rev. D 31, 725 (1985). 
[6] W. Israel; Ann. Phys. (N. Y.) 100, 310 (1976). 

[7] A. L. Garcia-Perciante, A. Sandoval- Villalbazo and L. S. Garcia-Colin; Gen. Rel. and Grav. 
(2009), online first. 

[8] C. Cercignani y G. Medeiros Kremer; The relativistic Boltzmann equation: theory and appli- 
cations, Birkhauser, Berlin (2002) 

[9] A. Sandoval- Villalbazo, A. L. Garcia-Perciante, L.S. Garcia-Colin, submitted to Physica A 
(2009), [Arxiv:gr-qc:0805.4237] 

[10] R. L. Liboff and R. C. Liboff; Kinetic theory: classical, quantum and relativistic descriptions, 
Springer, Berlin (2003). 

[II] E.W. Kolb and M.S. Turner The early universe Addison- Wesley, Reading Massachusetts(1990) 



6 



